Bohmian Mechanics with Complex Action: 
A New Trajectory-Based Formulation of Quantum Mechanics 
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In recent years there has been a resurgence of interest in Bohmian mechanics as a numerical tool 
because of its local dynamics, which suggest the possibility of significant computational advantages 
for the simulation of large quantum systems. However, closer inspection of the Bohmian formulation 
reveals that the nonlocality of quantum mechanics has not disappeared — it has simply been swept 
under the rug into the quantum force. In this paper we present a new formulation of Bohmian 
mechanics in which the quantum action, S, is taken to be complex. This leads to a single equation 
for complex 5*, and ultimately complex x and p but there is a reward for this complexification — 
a significantly higher degree of localization. The quantum force in the new approach vanishes for 
Gaussian wavepacket dynamics, and its effect on barrier tunneling processes is orders of magnitude 
lower than that of the classical force. We demonstrate tunneling probabilities that are in virtually 
perfect agreement with the exact quantum mechanics down to lO"'^ calculated from strictly localized 
quantum trajectories that do not communicate with their neighbors. The new formulation may have 
significant implications for fundamental quantum mechanics, ranging from the interpretation of non- 
locality to measures of quantum complexity. 
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Ever since the advent of Quantum Mechanics, there 
has been a quest for a trajectory-based formulation of 
quantum theory that is exact. In the 1950's, David 
Bohm, building on earlier work by Madelung 1] and de 
BroglieJ3|, developed an exact formulation of quantum 
mechanics in which trajectories evolve in the presence of 
the usual Newtonian force plus an additional quantum 
force j3i]. In recent years there has been a resurgence of 
interest in Bohmian mechanics (BM) as a numerical tool 
because of its local dynamics, which suggests the possibil- 
ity of significant computational advantages for the simu- 
lation of large quantum systems _4, _5, _6, _7, _8, _9, _10, Jd]. 
However, closer inspection of the Bohmian formulation 
reveals that the non-locality of quantum mechanics has 
not disappeared — it has simply been swept under the 
rug into the quantum force. Particularly disturbing is the 
fact that for simple cases such as Gaussian wave packet 
dynamics of the free particle or the harmonic oscillator, 
where classical-quantum correspondence should be per- 
fect, the quantum force is not only non- vanishing but is 
the same magnitude as the classical force 12j. 

In this paper we present a new formulation of BM in 
which the quantum phase, 5', is taken to be complex. This 
leads to a single equation for the complex phase, as op- 
posed to coupled equations for real phase and real am- 
plitude in the conventional BM. Complex phase leads to 
equations of motion for trajectories with complex x and 
p but there is a reward for this complexification — a 
significantly higher degree of localization than in conven- 
tional BM. We demonstrate tunneling probabilities that 
are in virtually perfect agreement with the exact quan- 
tum mechanics down to 10~^ calculated from strictly lo- 
calized quantum trajectories that do not communicate 
with their neighbors. There is a superficial similarity 



with some earlier workjlj, |lj| on a time-dependent ex- 
tension of WKB, but the present approach, which we 
call Bohmian mechanics with complex action (BOMCA) 
is formally exact, and not semiclassical. 

The starting point of conventional BM formulation (in 
1-dimension) is the insertion of the ansatz 



"ipix, t) = A{x, t) exp 



':S{x,t) 



(0.1) 



in the time dependent Schrodinger equation (TDSE), 
where A{x,t), S{x,t) are real functions representing the 
amplitude and phase respectively. Separating the result 
into its real and imaginary parts, two PDE's are obtained 
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where V{x) is the potential of the system. The first equa- 
tion is referred to as the quantum Hamilton- Jacobi (HJ) 
equation; it differs from the classical HJ equation (the 
LHS) by the addition of a "quantum potential" Q = 
^im4l^- Defining a velocity field v{x,t) = Sx{x,t)/m 
the classical HJ equation yields Newton's equation of 
motion; the same process for the quantum HJ equation 
yields equations of motion for "quantum trajectories". 
Eq. (|0.3() can be reformulated as a hydrodynamic-like con- 
tinuity equation for probability flow, hence eas. H().2l) and 
(10.3(1 are referred to as the hydrodynamic formulation of 
quantum mechanics. The solution of the quantum hydro- 
dynamic equations along the quantum trajectories con- 
stitutes the conventional BM formulation. 

The starting point of the BOMCA formulation is the 



insertion of the ansatz[l. 

ip{x, t) = exp 



^5(x,t) 



(0.4) 



in the TDSE, where we allow the phase to be complex. 
This yields a single newly defined quantum HJ equation 



St + 7; — S,j, + V — - — Sxx, 
2m 2ni 



(0.5) 



where Q = —-^S^x is the new quantum potential. Note 
that there is no expansion in powers of h in the derivation, 
hence eq. (|0.5|) is an exact formulation of the TDSE which 
to the best of our knowledge has not been explored in the 
literature. 

In the spirit of conventional BM our aim is to solve 
ea. (|0.5|) in the Lagrangian approach, that is, along quan- 
tum trajectories. A quantum trajectory is defined by 



dx 



= v{x,t); v{x,t) = —Sx{x,t). (0.6) 



Due to the definition of x as time dependent in eq. H0.6|l we 
write the solutions of this equation as x(i; xq) where xq is 
the starting point of the trajectory. Unlike conventional 
BM, the complex value of S yields quantum trajectories 
x{t; Xq) that evolve in the complex plane. As a conse- 
quence, the new formulation requires analytic continua- 
tion of the wavefunction and the phase to the complex 
plane. We consider only analytic potentials V{x) therefor 
v{x,t) is analytic in regions which do not contain nodes 
of ip{x, t). To obtain an equation of motion for v(x, t) we 
take the spatial derivative of eq. (|0.5() and apply eq. (|0.6|) 
to obtain 

ih 

■m{vt + Wx) - —Vxx = -Vx{x). (0.7) 

Identifying the expression in the round brackets as a La- 
grangian time derivative (^ = ■§i + v-§^) oiv, transforms 
ea. H0.7(l to a Newtonian- like equation of motion for the 
velocity 



dv[x{t;xQ),t] 
It 




(0.8) 



where we identify F^, Fq as the classical and the quantum 
force respectively. The non-locality of quantum mechan- 
ics is manifested in the appearance of Vxx in the quantum 
force. This term prevents the first equation in IjU.BI) and 
ea. (|0.8|) from being a closed set. 

As in conventional BM, the main difficulty lies in es- 
timating the quantum force. We tackle this problem by 
taking iterated spatial partial derivatives of ea. HQ.8|l . The 
result can be written after a short manipulation as 
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where go = and g^ = E"=i (")w(^')w(""^'+^) forn > 1. 
The superscripts denote the order of a partial spatial 
derivative. The set of eas. H0.9|l and the first equation 
in HO. 6(1 are now an infinite but closed set that describes 
a fully local complex quantum trajectory. If these equa- 
tions tend to as n —>■ 00, we may obtain a numeri- 
cal approximation by truncating the infinite set at some 
n = N, thus replacing ea. (|0.8() with a system of iV -I- 1 
coupled ODEs. Since each individual equation in 1(0. 9|l 
depends on the consecutive v("+^-' , the truncation is done 



by setting v 



{N+l) 
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for the ti(")'s are given by 
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where we applied the definition from ((0.6(1 together with 
S ~ —ihlnip that follows from eg. ((0.4(1 . xq G C is an 
initial position of an arbitrary single trajectory. A sim- 
ilar iterative differentiation process was applied in con- 
ventional BMplla, lD| yielding a more complicated set of 
coupled equations of amplitude and phase derivatives. 

The relation v — Sx/m, identifies the phase field with 
an action field of the quantum trajectories. The equation 
of motion for the action along a trajectory is similar to 
its classical counterpart with the addition of the quantum 
potential 

dS[x{t;xn),t] 1 3 t/ , «^ ^n ii\ 
— = St+ vSx = -mv -V + — Ui:.(0.11) 

Having x,v and Vx, the action along a trajectory is ob- 
tained simply by adding the integral of ea. ((0.11() to the 
initial value S{xq,0) = —ih\'a[ip{xo,0)]. Inserting the 
action in eg. 1(0.4(1 yields the value of the wavefunction 
iJj[x{t;xo),t] — exp {j-S[x{t; xo),t]} at position x{t;xQ) 
in the complex plane. 

To obtain the wavefunction on the real axis at time tj, 
in principle we need to propagate a set of initial positions 
{xo } such that {xj(tf]Xo-)} G R at a specified time 
tf. Tracing back the initial positions from final positions 
resembles the computationally expensive "root search" 
problem familiar from the semiclassical literature. How- 
ever, here a set of initial positions can be readily ob- 
tained. Suppose the initial positions are restricted to a 
region where -7^x{tf;xo) ^ and that the mapping of 
initial positions xq to final positions xitf^xo) is an ana- 
lytic function. Then the inverse mapping a; i-+ xq (t/ ; x) 
is also analytic. Consequently we can write 

r^ dxJtfx') 
xo{tf,xt)^xo{tf;Xa) + J g^; ' dx', (0.12) 

where we are free to choose the final positions Xa, xt and 
the integration contour. For simplicity, suppose we have 
found an initial condition xo such that Xa = x(tf;xo) G 
R. Varying Xb e R and choosing the integration contour 



to be the real interval between Xa and Xb, we obtain from 
eg. (|0. 12(1 a curve xo{tf; Xb) of initial conditions that reach 
the real axis at tf . This can be translated to a numerical 
scheme for generating initial positions that map to the 
vicinity of the real axis, by writing an iterative discrete 
equation based on ea. HU.12|l 
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where we define Sxq. = xq- — xoj__-^ and Sxj = 
Xj{tf]XQ.) — Xj^i{tf;xo-_-^). Ax' is a small step along 
the real axis. Since the final positions are not exactly on 
the real axis an interpolation process is used to extract 
the complex phase along the real axis from the action 
values at the set of final positions {xj(tf] xq)}. 

As a numerical example we consider the one- 
dimensional scattering of an initial Gaussian wavepacket 
V'(x,0) = (2a/7r)i/4e[-«(^-^=)'+iP<=(^-^<=)] from an 
Eckart potential V{x) — _D/cosh (/3x). We take Xc — 
-.7, a = SOtt, D = 40, /3 = 4.32 and to = 30 (aU units 
are atomic units). In Fig ^ we depict several complex 
quantum trajectories for the case of translational energy 
E = p1/2ra = 0, for iV = 1. Note that the complex val- 
ues of X and p allow the trajectories to "tunnel" through 
the barrier centered at x = 0. In FigI21 (a) we com- 
pare the exact wavefunction at t = .85, -E = 50 with the 
BOMCA results for truncation at TV = 1, ..., 4. Note that 
the transmitted part of the wavefunction is nearly con- 
verged for TV = 1, suggesting that BOMCA will be very 
efficient for calculating tunneling probabilities. In FigEl 
(b) we consider the case of extremely deep tunneling — 
E = — and focus on the transmitted part of the wave- 
function. The method converges uniformly, and even for 
truncation at iV = 1 the agreement with the exact results 
is excellent. It is interesting to contrast the equations for 
A'' = 1 with classical mechanics. For A^ = 1 there is no 
quantum force: the equations of motion are precisely the 
classical equations of motion 
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albeit for complex x and v. There is however a nonzero 
quantum potential that gives an additional term to the 
action integral (cf. eg. (|0. 11(1 ). where the term Vx fulfils 
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feg. (|0.9|) with n = A^ = 1). It is interesting to note that 
the defining equations for A^ = 1 (eqs. IO.14IU.l51 and 
eg. (|0. 11(1 ) have appeared previously in the literature, in 
the context of semiclassical methods [l^ [iSj ; we empha- 
size that here they emerge only as a convenient trunca- 
tion to an otherwise exact quantum formulation. 

The asymptotic tunneling probability T{E) is calcu- 
lated by integrating the absolute square of the wave- 
function for a; > at a sufficiently long time. In FigO 




Initial trajectory 

— Consequent trajectories 

o Initial position 

* Final position 



-0.45!- 



-0.5 



0.5 

Real(x) 



1 



1.5 



FIG. 1: Several complex quantum trajectories for the scat- 
tering of a Gaussian from an Eckart barrier centered around 
X = where E = p^^/m = 0, t = 1 and iV = 1. The 
trajectories were obtained through the numerical scheme de- 
scribed in the text fea. H().18|l ). Note that the trajectories 
initiate at Re{xo) — —.7 = Xc and reach Re[xf{xo;tf)] > 
0, Im[xf{xo;tf)] ~ 0. Hence, these trajectories "tunnel" 
through the barrier. 



(a) and (b) we compare the exact tunneling probabili- 
ties as a function of E with the results obtained from 
BOMCA and conventional BM. The exact results were 
computed by a split operator wavepacket propagation. 
The BOMCA results were calculated by propagating 50 
complex quantum trajectories. The conventional BM re- 
sults were calculated using the numerical formulation de- 
veloped by Lopreore and Wyatty. The BOMCA for- 
mulation allowed the exploration of tunneling over the 
whole energy range, while the conventional BM formu- 
lation proved unstable at low energies {E ^ 4). More- 
over, the BOMCA results are significantly more accurate 
than the BM results for all energies below the barrier 
height {E < D), even using just the classical equations 
of motion (A = 1). Note the improvement in the accu- 
racy of the BOMCA results as A" increases, suggesting 
convergence to the exact quantum result. The issue of 
convergence will be studied more fully in future work. 

In summary we presented BOMCA, a novel formu- 
lation of Bohmian mechanics. This formulation yields 
simpler equations than conventional Bohmian mechan- 
ics (at the expense of complex trajectories). Moreover, 
BOMCA allows a direct and simple derivation of local 
uncoupled trajectories that may be used to reconstruct 
the wavefunction. The tunneling probabilities obtained 
by BOMCA for the scattering process were in excellent 
agreement with the exact results even in the extremely 
deep tunneling regime. We showed that even classical 
equations of motion with a small number of (complex) 




0.2 0.4 0.6 0.1 

Real(x) 



1.2 



FIG. 2: Exact wavefunction vs. BOMCA reconstructed 
wavefunction for the scattering of a Gaussian from an Eckart 
barrier. Plot (a) corresponds to t = .85 and E — 50. Plot (b) 
focuses on the transmitted part of the wavefunction for the 
case of extremely deep tunneling — i = 1 and E — 0. Note 
the convergence to the exact wavefunction as A'' increases. 



trajectories are sufficient to obtain very accurate results 
provided that an extra, nonclassical term is added to the 
action integral. We wish to acknowledge Prof. Edriss S. 
Titi for several useful discussions. 
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FIG. 3: (a) Comparison between the tunneling probabilities 
obtained by BOMCA, conventional BM and the exact results. 
The inset shows an enlargement of the results for E « 4.17, 
the last point for which the BM formulation was stable; for 
_E < 4 we could not obtain stable results from the conventional 
BM formulation, (b) Log of the relative divergence from the 
exact results. 
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